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When a system thermalizes it loses all local memory of its initial conditions. This is a general 
feature of open systems and is well described by equilibrium statistical mechanics. Even within a 
closed (or reversible) quantum system, where unitary time evolution retains all information about its 
initial state, subsystems can still thermalize using the rest of the system as an effective heat bath. 

Exceptions to quantum thermalization have been predicted and observed, but typically require 
inherent symmetries mm or noninteracting particles in the presence of static disorder [3H6]. The 
prediction of many-body localization (MBL), in which disordered quantum systems can fail to 
thermalize in spite of strong interactions and high excitation energy, was therefore surprising and has 
attracted considerable theoretical attention HBno]. Here we experimentally generate MBL states 
by applying an Ising Hamiltonian with long-range interactions and programmably random disorder 
to ten spins initialized far from equilibrium. We observe the essential signatures of MBL: memory 
retention of the initial state, a Poissonian distribution of energy level spacings, and entanglement 
growth in the system at long times. Our platform can be scaled to higher numbers of spins, where 
detailed modeling of MBL becomes impossible due to the complexity of representing such entangled 
quantum states. Moreover, the high degree of control in our experiment may guide the use of MBL 
states as potential quantum memories in naturally disordered quantum systems m- 


It is exceedingly rare in nature for systems to local¬ 
ize, or retain local information about their initial con¬ 
ditions at long times. In an important counterexample, 
Anderson demonstrated that localization can arise due 
to the presence of disorder, which can destructively scat¬ 
ter propagating waves and prevent transport of energy 
or particles [3]. Although this interference effect can be 
applied to generic quantum systems, most experimental 
work has been restricted to the narrow parameter regime 
of low excitation energies and no interparticle interac¬ 
tions [IHH]. 

Whether such localization persists in the more general 
case of arbitrary excitation energy and non-zero inter¬ 
particle interactions was theoretically explored by Ander¬ 
son [3], and more recently by others [MH]- This MBL 
phase is predicted to emerge for a broad set of interac¬ 
tion ranges and disorder strengths, though the precise 
phase diagram is not well known [ 12 ] since equilibrium 
statistical mechanics breaks down in the MBL phase and 
numerical simulations are limited to ~ 20 particles mm- 
Very recent experiments searching for MBL have mea¬ 
sured constrained mass transport and the breakdown of 
ergodicity in disordered atomic systems with interactions 

[mn]. 

Here we report the direct observation of MBL in a long- 
range transverse field Ising model with programmable, 
random disorder. This is a non-integrable model that 
cannot be mapped to noninteracting particles (a neces¬ 
sary condition for MBL HI) and we can easily tune the 
disorder strength and interaction range over a parameter 


space that exhibits this phenomenon. Our experiment is 
effectively a closed quantum system over the timescales 
of interest, since the system localizes approximately ten 
times faster than the coupling rate to the outside envi¬ 
ronment. 

Each of the effective spin-1/2 particles is encoded in 
the ^Si /2 \F = 0, mp = 0) and |F = 1, mp = 0) hyperfine 
‘clock’ states of a ^’^^Yb+ ion, denoted j),)^ and It)^, re¬ 
spectively m- We confine a chain of 10 ions in a linear 
rf Paul trap and apply optical dipole forces to generate 
the effective spin-spin coupling [16] of a disordered Ising 
Hamiltonian: 

iJising = X + f E + E 

i<j i i 

where aj (7 = x, z) are the Pauli matrices acting on 
the spin, Jij is the coupling strength between spins 
i and j, B is a uniform effective transverse field, Di is 
a site-dependent disordered potential, and h = 1 (Meth¬ 
ods). After the chain evolves for some time, we collect 
the state-dependent fluorescence on an intensified charge- 
coupled device camera for site-resolved imaging. This, in 
addition to our ability to perform high fidelity rotations, 
allows measurement of the single-site magnetization (tr/) 
(7 = X, y, z) as well as arbitrary spin correlation functions 
along any direction. 

Jij is a tunable, long-range coupling that falls off ap¬ 
proximately algebraically as Jij (x Jmax/\i — j r m, 
where Jmax is typically 27r(0.5 kHz). Here we tune a 
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FIG. 1: An Interacting spin model with random disor¬ 
der (a) Directly measured elements of the spin-spin coupling 
matrix Jij (Eq. (1)). The long range interactions decay as 
. (b) A specific instance of the random disordered 
field with a schematic illustration of the long-range interac¬ 
tions and (c) the random values of the disordered field for all 
30 instances of disorder for several different disorder strengths 
and for each ion. (d) The level statistics calculated from the 
measured spin-spin coupling matrix (a) and applied disorders 
(c) are Poisson-distributed (black line is the expected level 
spacings for a Poisson distribution), as predicted for a MBL 
system. 


between 0.95 and 1.81, although for most of the data 
a « 1.13. Moreover, we directly measure the complete 
spin-spin coupling matrix (Fig. demonstrating the 
long-range interactions required to exhibit MBL. 

The site-specific programmable disorder term Di is 
sampled from a uniform random distribution with Di € 
[—IF, IF]. The disorder is generated by site-dependent 
laser-induced Stark shifts (Methods), which also allow 
for preparation of the system into any desired product 
state. To ensure we observe the general behavior of the 
disordered Hamiltonian, we average over 30 distinct ran¬ 
dom instances of disorder (Fig. B^-c), which leads to a 
sampling error that is smaller than the features of interest 
in the data. 

An important signature of the MBL phase is mani¬ 
fested in the spectral statistics of adjacent energy levels 
of the Hamiltonian. In a thermalizing phase, these en¬ 
ergy splittings follow random-matrix level statistics due 
to level repulsion. However, in the MBL phase, this 
level repulsion is greatly suppressed since eigenstates typ¬ 
ically differ by multiple spins flips. As a result, the level 
spacing between adjacent energy eigenvalues are Poisson- 
distributed Using our directly measured spin-spin 

couplings and applied realizations for the strongest ex¬ 


perimental disorder IF = SJmax and B = 4:Jmaxi we cal¬ 
culate the distribution of adjacent energy level splittings 
and find them to be Poisson-distributed, as expected for 
a MBL state (Fig. m 

Before searching for evidence of localization in the sys¬ 
tem’s time evolution, we first find parameters that cause 
the measured state to thermalize in the absence of disor¬ 
der. We increase the transverse field B and look for con¬ 
ditions that result in the single-site magnetization along 
two orthogonal directions approaching and remaining at 
their thermal equilibrium values (Methods). 

Figure shows the measured dynamics of (erf) for 
B = 4Jmax and Di = 0 with the spins initialized in the 
Neel ordered state, jtititititDz along the ^-direction. 
This configuration has an energy equivalent to an infinite 
temperature thermal state, since the expectation value of 
the Hamiltonian is zero. At long times, each expectation 
value erf approaches zero, losing memory of the initial or¬ 
dering. As the transverse field B is increased, the system 
appears to thermalize more quickly and the level statis¬ 
tics approach those of random matrices rather than Pois- 
sonians, as expected for a generic thermodynamic system 
(Methods). 

When B ^ J, the Hamiltonian is effectively an XY 
model HHHH] and conserves because Ising pro¬ 

cesses that flip spins along the large field are energetically 
forbidden. Thus, being in a spin configuration with half 
of the spins up and half of the spins down maximizes the 
accessible energy states. In addition, the Neel state is 
never an eigenstate, even for H ^ J and IF ^ J, since 
the uniform B field at each site still allows spin exchange 
in the z-basis. 

If a system is thermal, the Eigenstate Thermalization 
Hypothesis (ETH) provides a general framework where 
observables reach the value predicted by the microcanon- 
ical ensemble This allows us to calculate the 

expected thermal value of the reduced density matrix 
given the Hamiltonian and an initial state (Methods). 
To further establish that the system is thermalizing, 
we measure the reduced density matrix for each spin. 
Pi =Tr^j^iyp, without applied disorder and B = 4Jmax 
as shown in Eig. In our experiment, the spins are 
initially prepared in a product state with high fidelity. 
However at long times, the measured reduced density 
matrices show that each of the spins are very close to 
the zero magnetization mixed state, implying the system 
has locally thermalized. 

We apply the random disordered potential, Di ^ 0, 
and observe the emergence of MBL as we increase the 
strength of disorder. Since the many-body eigenstates 
in the MBL phase are not thermal, transport of energy 
and spins is suppressed, and ETH fails. Thus, observ¬ 
ables will not relax to their thermal values [9] and there 
will be memory of the initial conditions evident in the 
single-site magnetization. When starting in the Neel or¬ 
dered state. Fig. -f shows the time evolution of erf for 
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FIG. 2: Emergence of a Many-Body Localized State, (a) shows the time-evolved single-site magnetizations (erf) (different 
colors represent different ions) for the Hamiltonian in Eq. (1) and with B = 4:Jmax with no applied disorder {Di = 0). The 
initial-state reduced density matrices for ions 1 and 10 show the spins start in a product state along the 2 : direction. The 
time-averaged reduced density matrices for Jmaxt > 5 agree with the values predicted by the ETH, implying the system 
has thermalized locally, (b-e) As the disorder strength increases the spins retain more information about their initial state, 
indicating a transition towards MBL. (f) shows the dynamics of (nf) for the strongest applied disorder, W = 8Jmax- The initial 
and steady-state time-averaged reduced density matrices for ions 1 and 10 now show that information is preserved about the 
initial spin configuration at the end of the evolution. Statistical error bars (1 s.d.) are smaller than the data points. 


different disorder strengths. The frozen moments of the 
spins increase with increasing disorder as the emergent 
integrals of motion become more strongly localized uni. 

With the maximum applied disorder, W = iJmax, 
we measure the single-spin reduced density matrix for 
the initial state and the averaged matrix for Jmaxt > 5. 
In this case, localization of the spins leads to a marked 
difference in the measured and thermal reduced density 
matrices, indicating memory of the system’s initial con¬ 
ditions and a breakdown of ETH. 

To quantify the localization, we measure the normal¬ 
ized Hamming distance (HD) |53]: 

m = ^ ^ E (^oi w<(o) \^o) (2) 

i 

which gives the number of spin flips away from the initial 
state, normalized by the length of the chain. At long 
times, the HD approaches 0.5 for a thermalizing state 
and remains at 0 for a fully localized state. In Fig. ii, 
we measure that the long-time HD is 0.5 in the absence 
of disorder, and becomes smaller as the disorder strength 
is increased and the system more strongly localizes. 

Figure [S]] shows that for finite but weak disorder, the 
time-averaged HD for Jmaxt > 5 is unchanged, indicat¬ 
ing no localization. However, once the random field is 
sufficiently strong we observe a crossover from a ther¬ 
malizing to a localized state. Once in this regime, the 


system becomes more localized with increasing disorder 
strength. 

There is great theoretical interest in mapping the MBL 
phase diagram with respect to interaction range and dis¬ 
order strength [I1IMI21I- We have taken the first steps 
towards this goal by measuring a change in the time- 
averaged HD for W = SJmax and Jmaxt > 5 as we ad¬ 
justed the interaction range, 0.95 < a < 1.81 (Fig. §=)■ 
For shorter-range interactions, the system appears more 
localized, since the state approaches a fully-localized An¬ 
derson insulator as a —>■ 00 . 

Although there are predictions of a many-body delocal¬ 
ization transition at a = 1.5 for Hamiltonians similar to 
ours, we did not observe this effect as we tuned a across 
this boundary. The lack of a sharp transition, along with 
the presence of MBL states for a < 1.5, may be due to 
finite size effects (Methods). As this system is scaled to 
many dozens of spins, it will allow better study of the 
phase transition and mapping of the phase boundary in 
a regime where numerics are intractable. 

A hallmark of MBL is the characteristic growth of en¬ 
tanglement under coherent time evolution m , though its 
experimental observation has been elusive so far. In An¬ 
derson insulators without many-body interactions, the 
entanglement production from weakly entangled initial 
states shows a quick saturation after a sharp transient 
regime. However, in MBL systems a long-time growth 
sets in, which is logarithmically slow for short-range in- 
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FIG. 3: Hamming distance (HD), (a) The Hamming Dis¬ 
tance (HD) exhibits time dynamics that reach their steady- 
state values after Jmaxt ~ 5. For increasing disorder, the 
system becomes more strongly localized, and the steady-state 
Hamming Distance decreases, (different colors represent dif¬ 
ferent disorder strengths), (b) The steady-state HD with re¬ 
spect to the strength of the random potential indicates the 
state is not localized for small disorder, but after the random 
field is sufficiently strong it becomes more localized with in¬ 
creased disorder, (c) The system becomes less localized in the 
presence of longer-range interactions (smaller a). Error bars, 
1 s.d. 
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FIG. 4: Quantum Fisher Information (QFI). (a) The 

time evolution of the QFI for no disorder where there is no 
long-time growth of entanglement. The shaded area indicates 
the fast initial growth of QFI that follows a Lieb-Robinson- 
type bound, (b) The long-time logarithmic growth of the QFI 
for the applied disorder oi W — {6,8) Jmax is a lower bound 
for the entanglement growth in the MBL state. Black lines 
are logarithmic fits to the data. Statistical error bars (1 s.d.) 
are smaller than the data points. 


We have presented the experimental observation of 
MBL states in a quantum simulator with long-range in¬ 
teractions and random disorder. In a system whose level 
statistics predict MBL behavior, we observe a crossover 
between thermal and localized regimes as we increase 
the strength of applied random disorder, and we wit¬ 
ness a long-time growth of entanglement in the localized 
state. Our experimental platform is well suited for study¬ 
ing deep and intractable questions about thermalization 
and localization in quantum many-body systems. 
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METHODS 

Generating the effective Hamiltonian 


We generate spin-spin interactions by applying spin- 
dependent optical dipole forces to ions confined in a 3- 
layer linear Paul trap with a 4.8 MHz radial frequency. 
Two off-resonant laser beams with a wavevector differ¬ 
ence Ak along a principal axis of transverse motion glob¬ 
ally address the ions and drive stimulated Raman transi¬ 
tions. The two beams contain a pair of beatnote frequen¬ 
cies symmetrically detuned from the resonant transition 
at vq = 12.642819 GHz by a frequency /i, comparable to 
the transverse motional mode frequencies. In the Lamb- 
Dicke regime, this results in the Ising-type Hamiltonian 
in Eq. (1) [161131115^ with 


N 

m—1 




( 3 ) 


where H is the global Rabi frequency, uir = HAk'^ / {2M) 
is the recoil frequency, 6^ „j is the normal-mode matrix 
|33) . and ujm are the transverse mode frequencies. The 
coupling prohle may be approximated as a power-law de¬ 
cay Jij « Jm,ax/\i — j|“) where in principle a can be 
tuned between 0 and 3 by varying the laser detuning /j, 
or the trap frequencies uJm naEi]. In this work, a was 
tuned approximately between 0.95 and 1.81 by changing 
/r. By asymmetrically adjusting the laser beatnote de¬ 
tuning fj, about the carrier by a value of B we apply a 
global Stark shift that can be thought of as a uniform 
effective transverse magnetic field of {Bj2)af . 

We generate the effective disorder by applying a site- 
dependent Stark shift using a single 355nm laser beam 
that is focused down tightly to a 1/e^ waist of ~ l.S/rm. 
The ion separation is ^ 2.5^m, thus the crosstalk be¬ 
tween ions is negligible with a measured ratio of nearest- 
neighbor Rabi frequencies of ^ 20 : 1. We use an acousto 
optic modulator (AOM) with a full width at half maxi¬ 
mum bandwidth of «100 MHz to apply the Stark shift 
to each ion. The AOM is not imaged onto the ions, so 
that driving the AOM with different frequencies allows 
the position of the beam to be scanned over the length of 
a 10 ion chain, ^ 20/rm. The Stark shift is proportional 
to P. Thus, to achieve larger applied Stark shifts, we 
raster through the AOM drive frequencies corresponding 
to addressing each ion with a total cycle time of ~ 5fis in¬ 
stead of applying them simultaneously. Since we cannot 
control the sign of the site-specihc Stark shift, to center 
the disorder strength around the global transverse field, 
we adjust the asymmetric detuning by W Jmax/‘2- 
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FIG. 5: Measured thermalization in the transverse field Ising model. The upper panels show the time dynamics of 
(erf) (different colors represent different ions) for 10 spins prepared in the random product state |44.4't'U4't4't)a;) for different 
transverse magnetic field strengths. For B — 0 the spins are in a eigenstate and do not thermalize. However, as the strength of 
B is increased the system begins to thermalize more robustly and qnickly. The lower panel plots the time evolution of (erf) with 
10 spins prepared in the Neel ordered in the z direction for different transverse magnetic field strengths. We conclude that the 
system is in the thermalizing regime for B — iJmax since we observe thermalizing behavior along two orthogonal directions. 
Error bars are 1 standard deviation of statistical error. 


Measuring the spin-spin coupling matrix 

In order to observe the dynamics between just two of 
the ions in the chain, we shelve the other spins out of 
the interaction space. This is done by performing a tt 
rotation between ^Si /2 \F = 0,mp = 0), and one of 
the Zeeman states, ^Si /2 \F = l,mp = —1), while shift¬ 
ing the two ions of interest out of resonance by apply¬ 
ing a large Stark shift with the individual addressing 
beam. We then apply our Hamiltonian which now acts 
only on the two ions left in the interaction space and de¬ 
termine the elements of the spin-spin coupling matrix by 
fitting the measured interaction Rabi flopping frequency 
between each pair of spins. 


Arbitrary product state preparation 

State initialization starts with optically pumping the 
spins with high-fidelity to ■ ■ ■)z- Then we perform 
a global 1^12 rotation to bring the ions to ' ")x- 
this point we apply a Stark shift with the individual ad¬ 
dressing beam to the spins that are to be flipped and 
allow the chain to evolve until these ions are tt out of 
phase with rest of the ions. This, along with our ability 
to perform high fidelity global rotations, allows for the 
preparation of any arbitrary product state along any di¬ 
rection of the Bloch sphere. Individual spin flips can be 
achieved with a fidelity of ~ 0.97, while arbitrary state 
preparation can be done with a fidelity of ^ (0.97)-^, 


where N is the number of spins flipped with the individ¬ 
ual addressing beam. 


Determining a Set of Thermalizing Parameters 

Extended Data Figure shows the time evolution of 
(tjf) for different values of B for the spins initialized in 
the randomly chosen product state With¬ 

out a transverse field, the spins are in an eigenstate of 
the Ising interaction and undergo no evolution. Once a 
transverse field is added the individual spins begin to lose 
memory of their initial conditions and as its strength is 
increased, the ions thermalize faster and more robustly. 

To confirm the system is thermalizing, we measure 
the time evolution of the single site magnetization, (erf), 
along an orthogonal direction for different strengths of 
the transverse magnetic field starting with the spins ini¬ 
tialized in the Neel ordered state. As seen in Extended 
Data Fig. the spins have lost information about their 
initial conditions in the z direction for all values of B. 

We calculate the spectral statistics of adjacent energy 
levels for the Hamiltonian and find they are not Poisson 
distributed for B = ‘iJmax and Di — 0 indicating that 
with no applied disorder, the system is not in a localized 
phase. Furthermore, one can determine if a system is in a 
thermal or localized regime by finding the correlation be¬ 
tween adjacent energy splittings by calculating the ratio 
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FIG. 6: Thermalizing level statistics. The calculated 
value of (r) with respect to B saturates close to the predicted 
value for a random-matrix distribution (dashed black line) 
implying that the Hamiltonian is in the thermal phase for 
sufficiently large B. 


of two consecutive gaps [5]: 

_ min{5n,5n-i) 
max{dn,bn-i\ 

where 5n = En+i — > 0. For a localized phase, 

where one expects a Poisson energy spectrum, the prob¬ 
ability distribution of this order parameter is given by 
Pp{r) = 2/(1 -I- r)^ and thus (r) « 0.39. For energy level 
spacings following a random-matrix as predicted for a 
thermalizing regime, we calculate (r) « 0.53 for a chain 
of 10 spins. Extended Data Figure shows that (r) sat¬ 
urates to the expected value for a random matrix dis¬ 
tribution, indicating that the Hamiltonian is thermal for 
sufficiently large B. 
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FIG. 7: Comparison of the experimental data (green 
dots) with exact numerical simulations (thick blue 
lines) for QFI. The solid straight lines represent logarith¬ 
mic fits to the numerical (light blue) and experimental data 
(orange). Deviations from the ideal coherent dynamics due 
to decoherence in the experimental setup lead to a reduction 
of entanglement. Importantly, this suggests that the loss of 
purity does not generate a false positive for entanglement. 


Calculating the density matrix expected by the 
Eigenstate Thermalization Hypothesis 


Given a Hamiltonian and an initial state {tpo), the cor¬ 
responding energy is (V’ol ^ IV^o)- For a thermalizing sys¬ 
tem satisfying ETH this energy should be equal to the 
classical energy: 


Tr[i7e-'^^] 

Tr[e-/5^] 


(5) 


for the appropriate /3 = l/{kBT). When partitioning 
the entire system into subsystems A and B, with the size 
of A much smaller than B (perhaps even a single spin), 
then, the density matrix on site A at long times can be 
approximated by: 


PA Tr[e-/3^^] 


( 6 ) 


Since we start in the Neel ordered state, the initial 
energy given the Hamiltonian in Eq. (1) is equal to zero, 
(V'ol H I'i/o) = 0. Equating this to the right hand side of 
Eq. (§ and solving for /3 gives /3 = 0, or T = oo. Using 
this P in Eq. ([^ gives a value for any reduced thermal 
density matrix of: 


/ 1/2 0 \ 

1 0 1 / 2 ) 

in agreement with the measured reduced density matrices 
in Eig. 2a. 


Quantum Fisher Information 

The quantum Fisher information (QFI) has recently 
been shown to witness genuinely multipartite entangle¬ 
ment |34l 135] . From a quantum metrology perspective, 
the QFI quantifies the sensitivity of a given input state 










to a unitary transformation generated by the her- 
mitian operator O. In a pure state, it is given by [53] 

Fq = 4iAdf = 4i{6^) - (0)2). (7) 

For a local operator O = (where the difference 

between largest and smallest eigenvalue of Oi is 1 ), the 
QFI witnesses entanglement as soon as 

/q = Fq/7V>1. (8 ) 

To characterize the growth of entanglement out of the 
initial Neel state, the natural choice of the generator O is 
the staggered magnetization, O = 
markably, this QFI is proportional to the variance of the 
Hamming distance given by Eq. (2) of the main 
text, 

Fq = 4N^iAVf = , 

( 9 ) 

when associating 'D{t) = {V{t))^ with V = I/(2A^)[I — 

The QFI as defined in Eq. Q assumes a pure state, 
i.e., that time evolution is purely unitary. Eor mixed 
states, the QFI cannot be expressed as a simple expecta¬ 
tion value of the operator O [53]. In general, decoherence 
reduces the purity of the system’s state over experimental 
time scales. To show that the measured increase of Fq 
as defined in Eq. is indeed due to coherent dynamics, 
we compare to numerically exact calculations for a uni¬ 
tary time evolution using the experimental parameters. 
As can be seen in Eig. the experimental data is always 
below the theoretical prediction for a unitary time evo¬ 
lution. The loss of purity, therefore, does not generate 
a false positive indicator of entanglement in our system. 
Moreover, in the time window following the strong initial 
rise, the increase of Fq for both theory and experiment 
is consistent with a logarithmic growth (see below). 

To study how the localization behavior changes with 
system size, we performed a numerical finite-size scaling. 
In order to obtain a well-behaved scaling, we use the Kac 
prescription |33|, i.e., we adjust the couplings as Jij = 
where Af = (N - 1)~^ J2i<j I'i-- j\~°‘■ 
Note that using this prescription the fundamental energy 
scale J differs by about a factor of 2 from Jmax, the value 
used in the main text. 

Eor a > 2, the disordered long-range Ising Hamiltonian 
shows MBL behavior at large disorder m- In Fig. 1^ we 
plot the dynamics of the QFI for a = 3, where it can be 
seen that the system displays a characteristic long-time 
growth of entanglement. In particular, within a time win¬ 
dow 2“ < < J < 3“ where only next-nearest-neighbor in¬ 
teractions are relevant and next-to-next-nearest-neighbor 
interactions are still negligible, one expects the system to 
behave essentially as a nearest-neighbor Ising model with 
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FIG. 8: QFI from exact diagonalization {a = 3). 

When subject to disorder, the QFI of the staggered magneti¬ 
zation shows a characteristic growth of entanglement (bine 
lines; from dark to light: N = 4,6,..., 14 averaged over 
80000, 20000,..., 200 disorder realizations). This growth is 
absent in a theory of free fermions with long-range hopping 
and pairing (green dashed lines with A = 14 (dark green) 
averaged over 10000 realizations and N = 100 (light green) 
averaged over 1000 realizations). Inset: In a time window 
dominated by next-nearest neighbor interactions, 2“ < t J < 
3“, one observes a characteristic logarithmic entanglement 
growth, expected for a MBL system with short-range interac¬ 
tions. 


a weak next-to-nearest-neighbor coupling. For such a sys¬ 
tem, a logarithmic growth of entanglement is expected, 
as we indeed find in that regime, see inset in Fig.More¬ 
over, in Fig. we compare our numerical results for the 
spin chain to the appropriate long-range free fermionic 
theory (see below), which shows a quick system-size in¬ 
dependent saturation of entanglement without a further 
growth. Therefore, we conclude that the observed growth 
of the QFI is not possible in a quantum system without 
many-body interactions, thus giving a clear signature for 
true MBL behavior. 

The situation is more complex at a = 1.13. For the 
case H = 0, it has been predicted that within the range 
1 < a < 2 delocalized behavior could be expected in the 
thermodynamic limit m- As seen in Fig. for system 
sizes A < 12 the model displays all essential signatures of 
MBL, as found for 0 = 3. Performing finite-size scaling, 
however, indicates a crossover to potentially ergodic be¬ 
havior for larger chains. The growth of the QFI is cut off 
for larger N followed by a decrease of the entanglement 
on longer time scales. The time scale for the crossover be¬ 
comes smaller for increasing system size, which could in 
the thermodynamic limit potentially remove the charac¬ 
teristic entanglement growth completely. Hence, the sys¬ 
tem might become ergodic, but, unfortunately, the exact 
diagonalization quickly reaches its limits at this point. 
Here, scaling our quantum simulator to larger system 
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FIG. 9: QFI from exact diagonalization {a = 1.13 and 

W/J = 3). Same color coding as in Fig.[^ (a) For small sys¬ 
tem sizes N = 4,..., 10, we observe the characteristic growth 
of entanglement as for a = 3 (compare with Fig. [^. (b) For 
slightly larger systems {N > 10), a crossover appears at large 
times followed by a decrease of entanglement. This behavior 
might point towards a crossover to ergodicity in the thermo¬ 
dynamic limit. Importantly, however, for the experimentally 
relevant system size of = 10, we still find a time window 
consistent with a logarithmic growth of entanglement (see in¬ 
set). Notice the difference of about a factor of 2 in the time 
scale with respect to the main text. 

sizes could thus resolve a difficult open question, namely 
of the existence of ergodicity in the range 1 < a < 2. 
However, we would like to emphasize that the essential 
features of MBL phases in long-range disordered Ising 
chains are nevertheless captured by the considered ex¬ 


perimental system at N = 10. In particular, we still find 
a time window consistent with a logarithmic growth of 
entanglement, see inset in Fig. 

To show that the entanglement growth is truly due to 
interactions, we also compare the exact data to a close 
approximation of H, Eq. (1) of the main text, with a 
non-interacting theory. Using the Jordan-Wigner trans¬ 
formation, a~ —>■ with the phase of the string 

operator 9j = the Hamiltonian Eq. (1) can 

be mapped to a fermionic theory with annihilation and 
creation operators Cj and c| , respectively, 

i<j i 

( 10 ) 

If Jij contained only nearest-neighbor interactions, this 
Hamiltonian would become equivalent to a free fermionic 
theory. Eor general Jij, however, the string operators 
generate interactions between the fermions. Over short 
times, and especially in a localized regime, the phases 
Oj are dominated by their initial values, i.e., it is a 
good approximation to replace (for the initial Neel state) 
Oj —>■ TTX)j<i((~l)'^ + l)/2 in the Hamiltonian. This 
replacement amounts to approximating H hy a non¬ 
interacting fermionic theory with long-range hopping and 
pairing. The QFI for that case is included in Figs. 
and As one can see, the QFI quickly saturates to 
values below /q = 1. The experimentally and numer¬ 
ically observed further growth of the QFI is thus truly 
due to interactions, and cannot be captured within a free 
fermionic theory, even with long-range hopping. 












